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Classical and quantum annealing is discussed for a kinetically constrained chain of N non- 
interacting asymmetric double wells, represented by Ising spins in a longitudinal field h. It is 
shown that in certain cases, where the kinetic constraints may arise from infinitely high but vanish- 
ingly narrow barriers appearing in the relaxation path of the system, quantum annealing exploiting 
the quantum-mechanical penetration of sufficiently narrow barriers may be far more efficient than 
its thermal counterpart. We have used a semiclassical picture of scattering dynamics to do our 
simulation for the quantum system. 
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I. INTRODUCTION 

Here we demonstrate the effectiveness of quantum an- 
nealing [1] in the context of Kinetically Constrained Sys- 
tem (KCS) [2]. These KCSs are simple model systems 
which have trivial ground state structures and static 
properties, but a complex relaxation behaviour due to 
some explicit constraints introduced in the dynamics of 
the system [2]. Such systems are very important to un- 
derstand how much of the slow and complex relaxation 
behaviour of a glass can be attributed to its constrained 
dynamics alone, leaving aside any complexity of its en- 
ergy landscape structure. 

It has been demonstrated in certain models with en- 
ergy barriers [1,3] that one can effectively appoint quan- 
tum fluctuations (instead of thermal ones) to anneal a 
glassy system towards its ground state. In the method of 
quantum annealing, one introduces quantum fluctuations 
by including a term in the Hamiltonian due to tunnelling 
field, that does not commute with the original (classi- 
cal) Hamiltonian, and thus generate tunelling probabili- 
ties between the eigenstates (classical configurations) of 
the original (classical) Hamiltonian. The introduction 
of such a quantum tunnelling is supposed to make the 
infinitely high but infinitesimally narrow barriers trans- 
parent to the system. This allows transitions between 
different configurations classically trapped between such 
infinite barriers. In other words, it is expected that ap- 
plication of a quantum tunnelling term will make the free 
energy landscape ergodic, ie the system will consequently 
be able to visit any configuration with finite probability 
(Ray et al [1]). Finally, of course, the quantum tunnelling 
term is to be tuned to zero to get back the ground state 
of the classical Hamiltonian. 

To study quantum annealing [1,3] in a representative 
KCS it has to be appropriately generalised. The KCSs 
studied so far are all classical and the constraints are 
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absolutely unsurpassable. To incorporate the quantum 
effect we first visualize that these constraints originate 
from infinitely high energy barriers, so that the clas- 
sical system remains unable to jump over such a bar- 
rier at any finite temperature. Then in the quantum 
version of such a system we consider the possibility of 
tunnelling through such barriers quantum mechanically 
in certain cases when the barrier width approaches zero 
fast enough so that the barrier becomes integrable. We 
specifically define here a quantum version of a classical 
one-dimensional directed KCS, known as the East model 
[4], and study the quantum relaxational behaviour and 
consequent annealing (to the ground state in the classical 
limit). 

The classical East model is basically a one-dimensional 
chain of non-interacting Ising ('up-down') spins in a lon- 
gitudinal field h, say, in downward direction. The ground 
state of such a system is trivially given by all spins down. 
A kinetic constraint is introduced in the model by putting 
the restriction that the z-th spin cannot flip if the (i — 1)- 
th spin is down. Such a kinetic constraint essentially 
changes the topology of the configuration space, since the 
shortest path between any two configurations differing by 
one or more forbidden flips, is increased in a complicated 
manner owing to the blockage of the 'straight' path con- 
sisting of direct flips of the dissimilar spins. Further, 
the constraint becomes more limiting as more spins turn 
down, as happens in the late approach to equilibrium. 
As a result, the relaxation processes have to follow more 
complex and lengthier paths, giving rise to exponentially 
large timescale (~ e x l T , where T is the temperature) 
[4]- 



II. MODEL 

To introduce the possibility of quantum tunnelling 
through infinite but integrable barriers representing clas- 
sical constraints, we start with a chain of asymmetric 
double wells, each with a particle localized within it. 
When the barrier (or step) between the two wells is pen- 
etrable, then if we initially prepare a wave packet in one 
well with sufficiently high expectation value of kinetic en- 
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ergy T, it will eventually tunnel to the other well. If there 
is no dissipation, then such to and fro motion between the 
two wells will persist, following successive elastic reflec- 
tions from the infinitely high outer boundaries of the dou- 
ble well. In the corresponding Ising spin representation 
this tunnelling dynamics of a particle (wave packet) be- 
tween the two wells is represented by a quantum mechan- 
ical spin flip dynamics introduced in the model through 
the inclusion of a transverse field term in the Hamilto- 
nian. Clearly, the value of the transverse field T depends 
upon the height and width of the energy barrier between 
the wells and the kinetic energy T of the wave packets 
[5]. In such a representation, our model reduces to a 
chain of non-interacting Ising spins (double wells) in the 
presence of a downward field h (proportional to the well 
asymmetry). The spin flip dynamics (flipping probabil- 
ities) in this model will however be calculated directly 
from a semiclassical picture of the motion of a particle in 
a bounded double well. The kinetic constraint is intro- 
duced by assuming that the i-th spin faces an infinitely 
high energy barrier between its two states (up and down) , 
when the (i— l)-th spin is down. As in the classical model, 
this barrier is absent when the (i — l)-th spin is up (see 
Fig. 1). When the dynamics is classical, the barriers are 
impenetrable and the spin at i-th site has to wait until the 
(i — l)-th spin flips to the up state. In the quantum ver- 
sion of the model considered here, we allow for tunnelling 
through such (classically impenetrable) infinite barriers 
for the flip of the i-th spin even when the (i — l)-th spin 
is down. The tunnelling probabilities come from the fol- 
lowing semiclassical picture of scattering of a particle in 
a double well with infinitely remote outer boundaries (w 
— > oo in Fig. 1). If a particle is put in one of the wells of 
such a double well with some kinetic energy (actually the 
expectation value) F, then it will eventually be scattered 
by the separator (a barrier or step) between the two wells. 
In such a scattering, there is a finite probability P that 
the particle manages to go to the other well. We calcu- 
late P from the simple picture of scatterings of a particle 
by one dimensional potentials as prescribed below. The 
minimum of the energy of the Ising chain (eqivalent to 
the potential energy of the chain of the double wells) triv- 
ially corresponds to the state with all the spins down, i.e., 
aligned along the longitudinal field h (where all the par- 
ticles are in their respective lower wells). However, if one 
starts with a random configuration and kinetic energy 
r is not sufficient for tunnelling to the upper well, then 
the system, more or less, will exhibit the zero temper- 
ature (energy minimization) relaxation behaviour of the 
classical East model, and will extremely slowly approach 
the ground state (i.e., the minimum of the potential en- 
ergy). For sufficiently high F, the system occasionally 
tunnels through the infinite barriers corresponding to the 
constraints and thus can take up some of the relaxation 
paths forbidden classically. However, at any nonzero T, 
the ground state (lowest potential energy state) will be 
mixed with higher potential energy eigenstates. To reach 
the ground state, we start with a very large initial value 



of r and then reduce it following an exponential sched- 
ule given by T = T cxp (— £/tq). Here t denotes the time, 
and tq sets the effective time scale of annealing. At zero 
temperature the slow spin flip dynamics occurs only due 
to the tunnelling (kinetic energy) term F, and hence the 
system ceases to have any relaxational dynamics in the 
limit r — > 0. It may be mentioned here that in absence 
of any analytical expression for the tunnelling probabil- 
ity in asymmetric case of the type discussed here, (see 
e.g., [6]), we employ the asymmetric barrier tunnelling 
probabilities available [7]. 



III. SIMULATION AND RESULTS 

We have employed the quantum transmission (flipping) 
probablities (cf. [7]) from a very elementary scatter- 
ing picture which is qualitatively adequate, though not 
strictly valid for the asymmetric double well (shown in 
Fig. 1(b)) because of its bound states and finite w. Fol- 
lowing are the flipping probabilities (P) for the i-th spin 
in different possible situations used in our Monte Carlo 
simulation: 

I. If the (i — l)-th spin is up and the i-th spin is also 
up, then P = 1. 

II. If the (i — l)-th spin is up and the i-th spin is 
down, then (a) P = for T < 2h, and (b)P = 

mm{l, 4[r(r - 2h)] 1/2 /(VT + Vr - 2h) 2 }, for r > 2h. 

III. If the (i — l)-th spin is down and the i-th spin is up 
then P = m»n{l, 4[T(T + 2h)} 1/2 /((VT + y/T + 2h) 2 + 
9 2 )}- 

IV. If the (i — l)-th spin is down and the i-th spin 
is up then (a) P = for T < 2h, and (b)P = 

min{ i, 4 [r(r- 2h)] 1/2 /((VT+ Vr - 2hf+ g 2 )} for r > 

2h {h and T denoting the magnitudes only). 
Here g = x 2 a, \ and a being respectively the height and 
width of the barrier representing the kinetic constraint. 
The above expressions for P are actually the transmission 
coefficients in respective cases of one-dimensional scatter- 
ing across asymmetric barrier or step (according to the 
form of the potential encountered in passing from one 
well to the other, see e.g., [7]). Application of the above 
scattering picture, even for the double wells in Fig. lb 
(which our simulation is based on) as discussed before, is 
of course an approximation. 

In our simulation, we take N Ising spins (<7j = 
±1, i = 1, N) on a linear chain with periodic bound- 
ary condition. The initial spin configuration is taken to 
be random such that magnetization m — (l/A^J^CTi 
is practically negligible (mj s=s 0). We then start with 
a tunnelling field r and follow the zero temperature 
(semi-classical) Monte Carlo scheme as mentioned above, 
using the spin flip probabilities P's appropriate for the 
four cases I-IV. Each complete run over the entire lat- 
tice is taken as one time unit and as time progresses, 
r is decreased from its initial value Tq according to 
T = Foe - */ 1 "". The results are shown in Fig. 2. It 
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FIG. 1: Potential energy wells for the spin at site i, when (i — l)-th spin is (a) up and (b) down, with the external field h in 
the downward direction and barrier height x ~ * °°- For the classical East model, flipping to right well in (b) is impossible (at 
any finite T). In the quantum model considered here, although \ ~ * °°> the wel l width a — > in such a way that x 2ft (= <?) is 
finite, giving a finite tunnelling probability for going to the right well in (b). The limit w — » oo helps utilising the scattering 
picture employed here. 



shows that for g = 100 and To = 100 the system freezes 
before reaching the ground state (m/ = 1) for low values 
of tq; say for tq = 2000. For a somewhat greater value, 
e.g., tq = 5000, the system is completely annealed to the 
ground state within about 4 x 10 4 time steps. However, 
for a much greater tq, like tq — 20000, the system of 
course anneals completely but consumes more time un- 
necessarily. These generic features remain the same for 
other higher values of g. We have also studied the de- 
pendence of annealing behaviour with the parameter g, 
which is actually a measure of how impenetrable is the 
infinite barrier representing the kinetic constraint. Com- 
putations were carried out to locate, for a given value 
of g, the minimum value of tq for which the system 
just anneals upto m/ = 0.8 (complete annealing requires 
prohibitively longer computer time for this comparative 
study). We call this minimum value (rg) TO j n . A bisection 
scheme was used to locate (TQ) m i n for different values of 
g starting for the same initial configuration. The inset in 
Fig. 2 shows that (Tg) m i„ increases fairly sharply with 
g (an empirical analysis shows {TQ) min ~ ag 2 + b, where 
a and b are constants) for g < 1000. This quadratic 
variation with g might be due to the specific functional 
form of P we have used. However, for even higher values 
of g, the slope is expected to decrease, and finally in the 
asymptotic limit g — » oo, the relaxation behaviour should 
converge to that of one with an unsurpassable kinetic 
constraint (like the classical East model). This asymp- 
totic convergence could not however be explored, since 
the required computational time becomes prohibitively 
long as g is increased further. 

We have also studied thermal annealing of the same 
(classical East) model for T = 0, following an expo- 
nentially decreasing temperature schedule given by T — 
T exp (— t/rc); tc being the time constant for the ther- 
mal annealing schedule and T the initial temperature. 
Here, when (i — l)-th spin is down, the flipping probabil- 
ity for the i-ih spin (~ exp— (x/T)) is negligible since \ 



is very large. Otherwise, it flips with probability P = 1 if 
it were in the up state, and with Boltzmann probability 
P = exp (—h/T) if it were in the down state. In Fig. 3 
we compare the results for the same order of initial value 
and time constant for T and T (barrier height \ i s taken 
to be 1000 in both the cases while g was taken to be 
100 in the quantum annealing case, or equivalently the 
barrier width a is taken to be of the order of 10 -4 ). We 
observe that to achieve a similar degree of annealing (at- 
taining a certain final magnetization m/), starting from 
the same disordered configuration, one typically requires 
much smaller tq compared to tc', typically, tq ~ 10 3 x tq 
for equivalent annealing (for similar optimal values of fi- 
nal order m/ ~ 0.92). For annealing with final order 
m/ ~ 1, we find tc ~ 10 4 x tq. This comparison of 
course depends on the value of g used (for the barriers) 
as shown in the inset of Fig. 2. 



IV. SUMMARY 

We have discussed here the annealing of a kinetically 
constrained Ising spin chain of N spins, starting from a 
disordered state (with negligible initial magnetization), 
to its (external field induced) fully ordered ground state. 
At any finite temperature T (in the classical model) the 
system takes an exponentially long time to relax to the 
ordered state because of the kinetic constraints, which act 
like an infinite potential barrier, depending on the neigh- 
bouring spin configurations. Quantum mechanically, this 
infinite barrier is taken to be penetrable, ie with finite 
tunnelling probability, depending on the barrier height \ 
and width a (a — > faster than x~ 2 ). The introduced 
noise, required for the annealing, is reduced following an 
exponential schedule in both the cases: T = T§e~ t l TC , 
T = r e~ t / r '5, with Tq w r . For our simulation for the 
quantum case, we have taken the tunnelling probabilities 
P (for cases I-IV) and employed them in a semi-classical 
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FIG. 2: Quantum annealing (T = 0) for g = 100, To = 100 and h = 1 is shown for different values of tq, for a chain of 5xl0 4 
spins (m averaged over the same set C of 10 initial configurations for each tq). The horizontal (dashed) line indicates the average 
(over the same set C) value of m that could be reached from the initial configurations by simply minimizing the energy following 
the downhill principle. In the inset, variation of (TQ) min with g is shown (by the points) for one given configuration. The error 



in (rg) m i„ is typically less than 0.5 % in each case. The continuous line in the inset shows a fit (tq) t , 
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FIG. 3: Comparison between classical and quantum annealing for a chain of 5 x 10 4 spins (for the same initial disordered 
configuration with mi ~ 10 -3 ). We show the results for tq = 1.8 x 10 2 (for quantum) and tc = 10 6 (for classical) with h = 1; 
a lower tc would not produce substantial annealing. Starting from the same initial values To = To = 100, (and g = 100 in the 
quantum case) we observe that classical annealing requires about 10 7 steps, whereas quantum annealing takes about 10 4 steps 
for achieving the same final order m/ ~ 0.92. 



fashion for the one dimensional spin chain considered. 
We observe that for similar achievement in final order 
(to/ ~ 0.92 starting from rrii = 1CP 3 ), tc ~ 10 3 tq for 
iV = 5 x 10 4 . For even larger order (to/ ~ 1), quan- 
tum annealing works even better (tc ~ 10 3 tq, for the 
same value of N). These comparison are for g = 10 2 and 
X = 10 3 for the constraint barriers. 

In this picture, we considered the collective dynamics 
of a many particle system, where each one is confined 



in a (field) induced asymmetric double well potential for 
which we considered only the low lying two states (the 
wave packet localized in one well or the other), repre- 
senting the two states (up and down) of an Ising spin 
discussed above. The tunnelling of the wave packet from 
one well to the other was taken into account by employing 
a scattering picture and we used the tunnelling probabil- 
ities as the flip probabilities for the quantum Ising spins. 
As such, the reported simulation for the one dimensional 
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quantum East model is a semiclassical one. It may be 
noted however that, because of the absence of inter-spin 
interaction, the dimensionality actually plays no role in 
this model except for the fact that the kinetic constraints 
on any spin depend only on the left nearest neighbour 
(directedness in one dimension). Hence the semiclassical 
one dimensional simulation, instead of a proper quantum 
Monte Carlo simulation (equivalent to a higher dimen- 
sional classical one [5]), is quite appropriate here. 
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